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Abstract 

Existence of metastable quarks of new generation can be embedded into phe- 
nomenology of heterotic string together with new long range interaction, which only 
this new generation possesses. We discuss primordial quark production in the early 
Universe, their successive cosmological evolution and astrophysical effects, as well as 
possible production in present or future accelerators. In case of a charge symmetry 
of 4th generation quarks in Universe, they can be stored in neutral mesons, doubly 
positively charged baryons, while all the doubly negatively charged "baryons" are 
combined with He-4 into neutral nucleus-size atom-like states. The existence of all 
these anomalous stable particles may escape present experimental limits, being close 
to present and future experimental test. Due to the nuclear binding with He-4 pri- 
mordial lightest baryons of the 4th generation with charge +1 can also escape the 
experimental upper limits on anomalous isotopes of hydrogen, being compatible with 
upper limits on anomalous lithium. While 4th quark hadrons are rare, their presence 
may be nearly detectable in cosmic rays, muon and neutrino fluxes and cosmic electro- 
magnetic spectra. In case of charge asymmetry, a nontrivial solution for the problem 
of dark matter (DM) can be provided by excessive (meta)stable anti-up quarks of 4th 
generation, bound with He-4 in specific nuclear-interacting form of dark matter. Such 
candidate to DM is surprisingly close to Warm Dark Matter by its role in large scale 
structure formation. It catalyzes primordial heavy element production in Big Bang 
Nucleosynthesis and new types of nuclear transformations around us. 



1 Introduction 

The question about existence of new quarks and/or leptons is among the most important 
in the modern particle physics. Possibility of existence of new (meta)stable quarks which 
form new (meta)stable hadrons is of special interest. New stable hadrons can play the role 
of strongly interacting dark matter [1-3]. This question is believed to find solution in the 
framework of future Grand Unified Theory. A strong motivation for existence of new long- 
living hadrons comes from a possible solution [4] of the "doublet-triplet splitting" problem 
in supersymmetric GUT models. Phenomenology of string theory offers another motivation 
for new long lived hadrons. 

A natural extension of the Standard model can lead in the heterotic string phenomenol- 
ogy to the prediction of fourth generation of quarks and leptons [5, 6] with a stable 4th 
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neutrino [7-10]. The comparison between the rank of the unifying group Eq (r = 6) and 
the rank of the Standard model (r = 4) can imply the existence of new conserved charges. 
These charges can be related with (possibly strict) gauge symmetries. New strict gauge 
U(l) symmetry (similar to U(l) symmetry of electrodynamics) is excluded for known par- 
ticles but is possible, being ascribed to the fermions of 4th generation only. This provides 
theoretic motivation for a stability of the lightest fermion of 4th generation, assumed to 
be neutrino. Under the condition of existence of strictly conserved charge, associated to 
4th generation, the lightest 4th generation quark Q (either U or D) can decay only to 4th 
generation leptons owing to GUT-type interactions, what makes it sufficiently long living. 

Whatever physical reason was for a stability of new hypothetical particles, it extends 
potential for testing respective hypothesis due to its implications in cosmology. Especially 
rich in this sense is a hypothesis on (meta)stable quarks of new family. It defines the goal 
of current work. 

As we will show, in the case when 4th generation possesses strictly conserved U(l)- 
gauge charge (which will be called y-charge), 4th generation fermions are the source of new 
interaction of Coulomb type (which we'll call further y-interaction). It can be crucial for 
viability of model with equal amounts of 4th generation quarks and antiquarks in Universe. 
The case of cosmological excess of 4th generation antiquarks offers new form of dark matter 
with a very unusual properties. Owing to strict conservation of y-charge, this excess should 
be compensated by excess of 4th generation neutrinos. 

Recent analysis [11] of precision data on the Standard model parameters admits exis- 
tence of the 4th generation particles, satisfying direct experimental constraints which put 
lower limit 220 GeV for the mass of lightest quark [12]. 

If the lifetime of the lightest 4th generation quark exceeds the age of the Universe, 
primordial Q-quark (and Q-quark) hadrons should be present in the modern matter. If this 
lifetime is less than the age of the Universe, there should be no primordial 4th generation 
quarks, but they can be produced in cosmic ray interactions and be present in cosmic ray 
fluxes. The search for this quark is a challenge for the present and future accelerators. 

In the present work we will assume that up-quark of 4th generation (U) is lighter than its 
down-quark (D). The opposite assumption is found to be virtually excluded, if D is stable. 
The reason is that D-quarks might form stable hadrons with electric charges ±1 {(DDD)~ , 
(Du) + , (DDD) + ), which eventually form hydrogen-like atoms (hadron (DDD)~ is com- 
bined with A He ++ into +1 bound state), being strongly constrained in surrounding matter. 
It will become more clear from consideration of [/-quark case, presented below. 

The following hadron states containing (meta)stable [/-quarks (U-hadrons) are ex- 
pected to be (meta)stable and created in early Universe: "baryons" (Uud) + , (UUu) ++ , 
(UUU) ++ ; "antibaryons" (UUU)—, (UUu)~~, meson (Uuf. The absence in the Uni- 
verse of the states (Uud), (Uu) containing light antiquarks are suppressed because of 
baryon asymmetry. Stability of double and triple U bound states (UUu), (UUU) and 
(UU u), (UUU) is provided by the large chromo-Coulomb binding energy (oc oiq CD • uiq) 
[13, 14]. Formation of these states in particle interactions at accelerators and in cosmic rays 
is strongly suppressed, but they can form in early Universe and cosmological analysis of 
their relics can be of great importance for the search for 4th generation quarks. 

We analyze the mechanisms of production of metastable U (and U) hadrons in the early 
Universe, cosmic rays and accelerators and point the possible signatures of their existence. 
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We'll show that in case of charge symmetry of U-quarks in Universe, a few conditions 
play a crucial role for viability of the model. An electromagnetic binding of (UUU)~~ 
with 4 He ++ into neutral nucleus-size atom-like state (O-Helium ) should be accompa- 
nied by a nuclear fusion of (Uud) + and 4 He ++ into lithium-like isotope [ 4 He(Uud)] in 
early Universe. The realization of such a fusion requires a marginal supposition concerning 
respective cross section. Furthermore, assumption of ?7(l)-gauge nature of the charge, asso- 
ciated to U-quarks, is needed to avoid a problem of overproduction of anomalous isotopes 
by means of an y-annihilation of U-relics ([ 4 He(Uud)], (UUu), (UUU), *He{UUU), 
A He{UUu), (Uu)). Residual amount of U-hadrons with respect to baryons in this case is 
estimated to be less than 1CT 10 in Universe in toto and less than 1CP 20 at the Earth. 

A negative sign charge asymmetry of U-quarks in Universe can provide a nontrivial 
solution for dark matter (DM) problem. For strictly conserved charge such asymmetry in U 
implies corresponding asymmetry in leptons of 4th generation. In this case the most of U in 
Universe are contained in O-Helium states [ 4 He(UUU)] and minor part of them in mesons 
U u. On the other hand the set of direct and indirect effects of relic U-hadrons existence 
provides the test in cosmic ray and underground experiments which can be decisive for 
this hypothesis. The main observational effects for asymmetric case do not depend on the 
existence of y-interaction. 

The structure of this paper is as the following. Section 2 is devoted to the charge sym- 
metric case of U-quarks. Cosmological evolution of U-quarks in early Universe is con- 
sidered in subsection 2.1, while in subsection 2.2 the evolution and all possible effects of 
U-quarks existence in our Galaxy are discussed. The case of charge asymmetry of quarks of 
4th generation in Universe is considered in Section 3. Section 4 is devoted to the questions 
of the search for U-quarks at accelerators. We summarize the results of our present study, 
developing earlier investigations [15, 16], in Conclusion. 

2 Charge symmetric case of U-quarks 

2.1 Primordial [/-hadrons from Big Bang Universe 

Freezing out of U-quarks 

In the early Universe at temperatures highly above their masses fermions of 4th genera- 
tion were in thermodynamical equilibrium with relativistic plasma. When in the course of 
expansion the temperature T falls down below the mass of the lightest [/-quark, m, equi- 
librium concentration of quark-antiquark pairs of 4th generation is given by 



where 54 = 6 is the effective number of their spin and colour degrees of freedom. We use 
the units h = c = k = 1 throughout this paper. 

The expansion rate of the Universe at RD-stage is given by the expression 




(1) 




(2) 
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where temperature dependence follows from the expression for critical density of the Uni- 
verse 

Pcrit ~ sttG ~ 9tot m ■ 

When it starts to exceed the rate of quark-antiquark annihilation 

Rann = ™4 (<Tv) , (3) 

in the period, corresponding to T = Tt < m, quarks of 4th generation freeze out, so that 
their concentration does not follow the equilibrium distribution Eq.(Q]) at T < Tf. For a 
convenience we introduce the variable 

T4 = —, (4) 

s 

where 

s = 2tt |p r 3 ^ 1MgtotaTh w 1.80 gZfsV^ns (5) 

is the entropy density of all matter. In Eq.© s was expressed through the thermal photon 
number density ra 7 = ^>r- T 3 and also through the baryon number density hb, for which 
at the modern epoch we have n B lod /n™ od = n « 6 • 10~ 10 . 

Under the condition of entropy conservation in the Universe, the number density of the 
frozen out particles can be simply found for any epoch through the corresponding thermal 
photon number density n T Factors g tot and g to ts take into account the contribution of all 
particle species and are defined as 

TA 4 7 ^ /Tj x4 



9tot = Y 9i ij;) + g Y 9 

i=bosons i=fermions 



T 



gtots = Y 9i[^) +\ Y 9i(Wl ■■ 



and 

k t J 8 . y \t ) 

i=bosons i=jermions 

where g± and Tj are the number of spin degrees of freedom and temperature of ultrarela- 
tivistic bosons or fermions. For epoch T <C m e « 0.5 MeV it is assumed that only photons 
and neutrinos with T v = (4/ll) 1 / 3 T give perceptible contribution into energy (until the 
end of RD-stage) and entropy (until now) densities so one has 

gZ d s « 3.91 gZf « 3.36. (6) 

For modern entropy density we have s mot i rj 2890 cm~ 3 . 

From the equality of the expressions Eq.© and Eq.© one gets 

I/O 

m/Tf 42 + lrL(g tot ' m p m (av)) 

a 2 

with m p being the proton mass and obtains, taking (av) ~ Q % D and gtot(Tf) = 
9tots(T f )=g f « 80-90, 

Tf « m/30 



4 



and 



r 4 = 




2.5 • 10 



-14 



m 



(7) 



s f W v ) g) l2 mpiT } (av) 



250 GeV 



Index "f ' means everywhere that the corresponding quantity is taken at T = Tf. Note, 
that the result Eq.fTJ), obtained in approximation of "instantaneous" freezing out, coincides 
with more accurate one if (av) and gj can be considered (as in given case) to be constant. 
Also it is worth to emphasize, that given estimation for relates to only 4th quark or 4th 
antiquark abundances, assumed in this part to be equal to each other. 

Note that if Tf > A = mo — tji, where mo is the mass of D-quark (assumed to be 
heavier, than [/-quark) the frozen out concentration of 4th generation quarks represent at 
Tf > T > A a mixture of nearly equal amounts of UU and DD pairs. 

At T < A the equilibrium ratio 



is supported by weak interaction, provided that /^-transitions (U — > D) and (D — > U) are in 
equilibrium. The lifetime of D-quarks, r, is also determined by the rate of weak (D — » U) 
transition, and at t 3> r all the frozen out DD pairs should decay to UU pairs. 

At the temperature Tf annihilation of U-quarks to gluons and to pairs of light quarks 
UU — ► gg, qq terminates and UU pairs are frozen out. The frozen out concentration is 
given by Eq.©. Even this value of primordial concentration of [/-quarks with the mass 
m = 250 GeV would lead to the contribution into the modern density 1mr±s mo d, which 
is by an order of magnitude less than the baryonic density, so that in the charge symmetric 
case [/-quarks can not play a significant dynamical role in the modern Universe. 

The actual value of primordial U -particle concentration should be much smaller due 
to QCD, hadronic and radiative recombination, which reduce the abundance of frozen out 
U -particles. y-Interaction can play essential role in successive evolution to be considered. 
It accounts for radiative recombination and plays crucial role in galactic evolution of U- 
hadrons. So, it will be included into further consideration which will be carried out for both 
sub-cases (with and without y-interaction). 

QCD recombination 



where a = 0.23 accounts for joint effect of Coulomb-like attraction due to QCD and y- 
interactions, formation of bound (UU) states is possible, in which frozen out Heavy quarks 
and antiquarks can annihilate. Effect of y-interaction is not essential here. 

Note that at T < I\ rate of (UU) annihilation in bound systems exceeds the rate of 
"ionization" of these systems by quark gluon plasma. So the rate of QCD recombination, 
given by [14, 15] 




At 



T <h = ma 2 /4 = 3.2 GeV 



m 



250 GeV 




(8) 



is the rate, with which abundance of frozen out U -quarks decreases. 
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The decrease of U -hadron abundance owing to UU recombination is governed by the 
equation 

= -3Fn 4 - n\ ■ (av) . (9) 

Using notation Eq.© and relation 

dT 

- dt =Hf> ( 10 ) 
which follows from Eq.© and is true as long as g to t ~ const, Eq.® is reduced to 

dr 4 = rl-sHT(av)dT, (11) 



where sHT = ^firg/Ab with g = gf ots /gtot = (for T > m e ) = g to ts = gtot- 

At To = I\ > T > Tqcd = T\, assuming in this period g = const = gj ps 17, the 
solution of Eq.dTTI) is given by 

rn f m \ 1 / 2 m 
r 4 = , ps 0.16 — ps (12) 

1 + ro VW m Pi It° M ^ t V 1 7 " mpi 

«1.6-10- 16 m 



250 GeV ' 

It turns to be independent on the frozen out concentration ro given by Eq.©. 

At T < Ijju < ma 2 /4 = 1.6 GeV 2 5o > GeV ' wnere effective constant a = CFa s — a y ~ 
(4/3) • 0.144 — 1/30 = 0.16 accounts for repulsion of the same sign y-charges, reactions 
U + U -> (UU) + g and U + (UU) (UUU) + g can lead to formation (£7[/)-diquark 
and colorless (UUU) "hadron" (as well as similar U bound states) in quark gluon plasma 
[13, 14]. However, disruption of these systems by gluons in inverse reactions prevents their 
effective formation at T > Iuu I '30 [14]. Therefore, such systems of U quarks with mass 
m < 700 GeV are not formed before QCD phase transition. 

Hadronic recombination 
After QCD phase transition at T = Tqcd ~ 150MeV quarks of 4th generation combine 
with light quarks into U -hadrons. In baryon asymmetrical Universe only excessive valence 
quarks should enter such hadrons. Multiple U states formation can start only in processes 
of hadronic recombination for U-quark mass m < 700 GeV what is discussed below. 

As it was revealed in [5, 6] in the collisions of such mesons and baryons recombination 
of U and U into unstable (UU) "charmonium -like" state can take place, thus successively 
reducing the [/-hadron abundance. Hadronic recombination should take place even in the 
absence of long range y-interaction of [/-particles. So, we give first the result without the 
account of radiative recombination induced by this interaction. 

There is a large uncertainties in the estimation of hadronic recombination rate. The 
maximal estimation for the reaction rate of recombination (av) is given by 

( av )^ J_ « 6 -10- 16 — (13) 

m~ s 

or by 

(av) ~i-«2-10- 17 — . (14) 

m~ s 
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The minimal realistic estimation gives [15] 



(av) « 0.4 • (r e// m 3 )- x /2(3 + log (T QCD /T eff )), (15) 

where T e jf = max {T, aym^}- 

Solution of Eq.dTTT) for (av) from the Eq.(fT3T> is given by 
Case A 

r 4 = =JL= « 1.0 • 10- 20 (16) 

14-r-n / n 9QCD mpi J-QCD 
x i ' Y 45 i% t% 

and it is (^p) 2 ~ 30 times larger for (av) from the Eq.(fl~4l): 
Case B 

r 4 = ^12= = «3.0-10- 19 . (17) 

1 _i_ ™ / n 9QCD mpi J-QCD 
1 + r y 45 mp mp 

For the minimal estimation of recombination rate ([TBI ) the solution of Eq. (fTTb has the form 

n = ^= ^= (18) 



45 m 

where in all the cases vq is given by Eq.(fT2l) and qqcd ~ 15. We neglect in our estimations 
possible effects of recombination in the intermediate period, when QCD phase transition 
proceeds. 

The solutions ([TBI ) and (TTTT ) are independent on the actual initial value of r 4 = ro, if 
before QCD phase transition it was of the order of (fl2l) . For the minimal estimation of the 
recombination rate (fT5l) the result of hadronic recombination reads 

CaseC 

-16 ( m \ 3//2 



r ^ 1 - 2 ' 10 UoGivJ • <19) 

As we mentioned above, for the smallest allowed mass of ?7-quark, diquarks (UU), 
(UU) and the triple U (and U) states (UUU), (UUU) can not form before QCD phase tran- 
sition. Therefore U-baryonic states (UUu), {UUU) and their antiparticles should originate 
from single U (and U) hadron collisions. The rate of their creation shares the same theo- 
retical uncertainty as in the case of (UU) formation, considered above. Moreover, while 
baryon (UUu) can be formed e.g. in reaction (Uud) + (Uud) — ► (UUu) + n, having no 
energetic threshold, formation of antibaryon (UUu) may be suppressed at smallest values 
of m by the threshold of nucleon production in reaction (Uu) + (Uu) — > (UU u) + p + tt + , 
which can even exceed UU binding energy. In further consideration we will not specify 
[7-hadronic content, assuming that (UUU), (UUu) and (Uu) can be present with appre- 
ciable fraction, while the content of residual £/-hadrons is likely to be realized with multiple 
U-states and with suppressed fraction of single U-states. Nevertheless we can not ignore 
single U-baryonic states (Uud) + because only reliable inference on their strong suppres- 
sion would avoid opposing to strong constraint on +1 heavy particles abundance which will 
be considered below. 

Radiative recombination 
Radiative UU recombination is induced by "Coulomb-like" attraction of U and U due to 
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their y-interaction. It can be described in the analogy to the process of free monopole- 
antimonopole annihilation considered in [17]. Potential energy of Coulomb-like interaction 
between U and U exceeds their thermal energy T at the distance 

do = J,- 

In the case of y-interaction its running constant a = a y ~ 1/30 [5]. For a -C 1, on the 
contrary to the case of monopoles [17] with g 2 /4it 3> 1, the mean free path of multiple 
scattering in plasma is given by 



A = (nay 1 ~ T 3 



2 \ _± 



Tm j a 3 T 



being A 3> do for all T < m. So the diffusion approximation [17] is not valid for our case. 
Therefore radiative capture of free U and U particles should be considered. According 
to [17], following the classical solution of energy loss due to radiation, converting infinite 
motion to finite, free U and U particles form bound systems at the impact parameter 

a w (T/m) 3 / 10 • d - (20) 

The rate of such binding is then given by 

lav) = 7ra 2 v « vr • (m/T) 9 / 10 • f-t « 

\m J 



(21) 



6 • 10~ 13 



a \ 2 / 300K \ 9 / 10 ^SOGeVY 1 / 10 cm 3 
l730j \ 



m 



The successive evolution of this highly excited atom-like bound system is determined 
by the loss of angular momentum owing to y-radiation. The time scale for the fall on the 
center in this bound system, resulting in UU recombination was estimated according to 
classical formula in [18] 

a 3 fm\ 2 a /m\ 21 / 10 1 

•- (22) 
m 



64tt ' V a J 64vr ' V T J 



V T J V250GeV; 

As it is easily seen from Eq.(l30l) this time scale of UU recombination r <C m/T 2 <C 
mpi/T 2 turns to be much less than the cosmological time at which the bound system was 
formed. 

The above classical description assumes a = m3 /i" T7 /i 3> ^ and is valid at T <C 
ma 20 / 7 [14]. 

Kinetic equation for U-particle abundance with the account of radiative capture on RD 
stage is given by Eq.dTTI). 
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At T < T rr = af /7 m « 10 MeV(m/250 GeV) (i^of° /7 the solution for the effect 



Vl/30 J 

of radiative recombination is given by 



r 4 ~ „ I 1 777777 ~ r (23) 



l + ro^Z^^) 1710 



with r taken at T = T rr equal to r 4 from Eqs.(fl6l),(fT7l or dT9b . 

Owing to more rapid cosmological expansion radiative capture of U -hadrons in expand- 
ing matter on MD stage is less effective, than on RD stage. So the result r 4 ro holds on 
MD stage with even better precision, than on RD stage. Therefore radiative capture does 
not change the estimation of [/-hadron pregalactic abundance, given by Eqs.(fT6l).(fT7l) or 

CE9J. 

On the galactic stage in the most of astrophysical bodies temperature is much less than 
T rr and radiative recombination plays dominant role in the decrease of U -hadron abundance 
inside dense matter bodies. 

U-hadrons during Big Bang Nucleosynthesis and thereafter 

One reminds that to the beginning of Big Bang Nucleosynthesis (BBN) there can be 
(Uud)+, (UUu) ++ , (UUU) ++ , (Uu)°, (UtJU)—, (UUu)~ states in plasma. We do 
not specify here possible fractions of each of the U-hadron species (i) in U-hadronic matter, 
assuming that any of them can be appreciable (r« < r 4 ). 

After BBN proceeded, the states (UUU) , (UUu) are combined with 4 He ++ due 
to electromagnetic interaction. The binding energy of the ground state can be estimated 
with reasonable accuracy following Bohr formulas (for point-like particles) 

I b = iZAZx f 2mA *1.5MeV, (24) 

where Zx = 2, Za = 2 and uia ~ 3.7 GeV are the charges of U-hadron and Helium and 
the mass of the latter. Cross section of this recombination is estimated as [19] 

. , 2 8 iry/2^a 3 Z 4 A Z 2 3.06 • 10~ 4 

(av) = j— f w (25) 

3 exp(4)mAv rriAT mAvrriAT 

Evolution of abundance of U-hadrons combining with He is described by equation 

dn (UUU) Q „ , . 

— = -3Hn {0m) - (av) n {t j W) n H e- (26) 

The term corresponding to disintegration of [(UUU)He] is neglected, since the energy of 
thermal photons is insufficient to disintegrate [(UUU) He] (the same for [(UUu) He]) in 
the ground state in this period. Following procedure Eqs.(|9lfTTI). we get 

r (uuu) = r (uuu)o ex P (~\/|f mpi f Q rHe ( av "l dT "j ~' (27) 

~ r (uuu)o ex P (-°- 6 • 1q12 
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where r# e = nue/s = Y p /4 ■ rj ■ n™ od j s mo d ~ 5.2 • 10~ 12 , g follows from Eq.© and 
To = 100 keV was taken. As one can see, Eq.(l27T) gives in this case strong exponential sup- 
pression of free (UUU) (the same for (UUu)), while neutral [(UUU) He] and [(UUu) He] 
states, being one of the forms of O-helium [16,20-26], catalyze additional annihilation of 
free [7-baryons and formation of primordial heavy elements [27]. New type of nuclear re- 
actions, catalyzed by O-helium, seem to change qualitatively the results of BBN, however 
(see Sec. [3] and arguments in [16, 20-27]) it does not lead to immediate contradiction with 
the observations. 

On the base of existing results of investigation of hyper-nuclei [28], one can expect that 
the isoscalar state = (U ud) + can form stable bound state with A He due to nuclear inter- 
action. The change of abundance of U-hyperons A^ owing to their nuclear fusion with 4 He 
is described by Eq. (126l27l) . substituting (UUU) <-> A^. Disintegration of [Ay He] is also 
negligible, since the period, when BBN is finished, is of interest (T < To -C I ([Ay He])). 

Cross section for nuclear reaction of question can be represented in conventional pa- 
rameterization through the so called astrophysical S -factor 

a = — exp(^ j, (28) 

where E = fiv 2 /2 with \i being reduced mass of interacting particles and v being their 
relative velocity. The exponent in Eq.d28l) expresses penetration factor, suppressing cross 
section, which reflects repulsive character of Coulomb force contrary to the case of (UUU). 
S-factor itself is unknown, being supposed S(E — > 0) — > const. Averaging av over 
Maxwell velocity distribution gives, using saddle point method, 

Av Q ■ S(E(v )) ( 
(av) -= exp — 

\ / ^ T 2T 

where v = ( 2 ™ z f* T f 3 . 

Calculation gives that suppression of free A^ on more than 20 orders of magnitude 
is reached at S(E) > 2 MeV • barn. S-factor for reaction of 4 He production is typically 
distinguished by high magnitudes from those of other reactions and lies around 5 — 30 MeV- 
barn [29]. However, reactions with 7 in final state, which is assumed in our case (Au + 
4 He — > [AjjHe] + 7), have as a rule S-factor in 10 4 times smaller. Special conditions 
should be demanded from unknown for sure physics of A;y-nucleus interaction to provide 
a large suppression of Ajj abundance. Such suppression is needed, as we will see below, to 
avoid contradiction with data on anomalous hydrogen abundance in terrestrial matter. The 
experimental constraints on anomalous lithium are less restrictive and can be satisfied in 
this case. 

y-plasma 

The existence of new massless U(l) gauge boson (y-photon) implies the presence of pri- 
mordial thermal y-photon background in the Universe. Such background should be in equi- 
librium with ordinary plasma and radiation until the lightest particle bearing y-charge (4th 
neutrino) freezes out. For the accepted value of 4th neutrino mass (> 50GeV) 4th neu- 
trino freezing out and correspondingly decoupling of y-photons takes place before the QCD 
phase transition, when the total number of effective degrees of freedom is sufficiently large 
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to suppress the effects of y-photon background in the period of Big Bang nucleosynthesis. 
This background does not interact with nucleons and does not influence the BBN reactions 
rate (its possible effect in formation and role of [ 4 He(UUU)] "atom" is discussed in [15]), 
while the suppression of y-photon energy density leads to insignificant effect in the speed- 
ing up cosmological expansion rate in the BBN period. In the framework of the present 
consideration the existence of primordial y-photons does not play any significant role in the 
successive evolution of [/-hadrons. 

Inclusion of stable y-charged 4th neutrinos strongly complicate the picture. Condition 
of cancellation of axial anomalies requires relationship between the values of y-charges of 
4th generation leptons (N, E) and quarks (U, D) as the following 

£yN = %E = —e y u/3 = —e y D/3. 

In course of cosmological combined evolution of U and N and y, "y-molecules" of kind U- 
U-U-N, where different U-quarks can belong to different U-hadrons (possibly bound with 
nucleus) should form. Such y-neutral molecules can avoid effect of U-hadrons suppression 
in the terrestrial matter, relevant in charge symmetric case, and lead to contradiction with 
observations, analysis of which is started now. UUU-N-type states will be considered in 
section [3] devoted to the charge-asymmetric case. 

2.2 Evolution and manifestations of U -hadrons at galactic stage 

In the period of recombination of nuclei with electrons the positively charged U -baryons 
recombine with electrons to form atoms of anomalous isotopes. The substantial (up to 
10 orders of magnitude) excess of electron number density over the number density of 
primordial [/-baryons makes virtually all [/-baryons to form atoms. The cosmological 
abundance of free charged [/-baryons is to be exponentially small after recombination. 

Hadrons (UUu), (UUU) form atoms of anomalous He at T ~ 2eV together with re- 
combination of ordinary helium. The states [(UUU)He], [(UUu)He], (Uu) escape recom- 
bination with electrons because of their neutrality; hadrons (Uud), if they are not involved 
into chain of nuclear transitions, form atoms of anomalous hydrogen. 

The formed atoms, having atomic cross sections of interaction with matter follow bary- 
onic matter in formation of astrophysical objects like gas clouds, stars and planets, when 
galaxies are formed. 

On the contrary, O-helium and (Uu) mesons, having nuclear and hadronic cross sec- 
tions, respectively, can decouple from plasma and radiation at T ~ 1 keV and behave in 
Galaxy as collisionless gas. In charge asymmetric case, considered in the next Section 
or in charge symmetric case without y-interaction O-helium and (Uu) mesons behave on 
this reason as collisionless gas of dark matter particles. On that reasons one can expect 
suppression of their concentration in baryonic matter. 

However, in charge symmetric case with y-interaction, the existence of Coulomb-like 
y-attraction will make them to obey the condition of neutrality in respect to the y-charge. 
Therefore owing to neutrality condition the number densities of U- and U -hadrons in as- 
trophysical bodies should be equal. It leads to effects in matter bodies, considered in this 
subsection. 
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U-hadrons in galactic matter 

In the astrophysical body with atomic number density n a the initial U -hadron abundance 
n uo = faO ' n a can decrease with time due to UU recombination. Here and in estimations 
thereafter we will refer to U-quark abundance as U-hadron one (as if all U-hadrons were 
composed of single U-quarks), if it is not specified otherwise. 
Under the neutrality condition 

njj = rifj 

the relative U -hadron abundance f a0 = njj jn a = rifj/n a is governed by the equation 

d ^ = -fl.n a .{av). (30) 
Here (av) is defined by Eq.(l2TT). The solution of this equation is given by 

f° = m fa ° i w - (31) 

1 + Jao ■ n a ■ (av) ■ t 

If 

n a -(av) -i> — , (32) 

JaO 



the solution (I3TI ) takes the form 



/« = ^XT- (33) 

n a ■ (av) ■ t 



and, being independent on the initial value, ?7-hadron abundance decreases inversely pro- 
portional to time. 

By definition f a o = fo/A a tom, where A a tom is the averaged atomic weight of the 
considered matter and /o is the initial U -hadron to baryon ratio. In the pregalactic matter 
this ratio is determined by from A) Eq.([T6l). B) Eq.(IT7b and C) Eq.([T9l) and is equal to 

10~ 10 for the case A, 
r -± = \ 3 • 10~ 9 for the case B, (34) 
1.2 • 10" 6 for the case C. 



n 



Here r\, « 10~ 10 is baryon to entropy ratio. 

Taking for averaged atomic number density in the Earth n a 10 23 cm -3 , one finds that 
during the age of the Solar system primordial U -hadron abundance in the terrestrial matter 
should have been reduced down to f a 10~ 28 . One should expect similar reduction of 
U -hadron concentration in Sun and all the other old sufficiently dense astrophysical bodies. 
Therefore in our own body we might contain just one of such heavy hadrons. However, as 
shown later on, the persistent pollution from the galactic gas nevertheless may increase this 
relic number density to much larger value (f a « 10 -23 ). 

The principal possibility of strong reduction in dense bodies for primordial abundance 
of exotic charge symmetric particles due to their recombination in unstable charmonium 
like systems was first revealed in [30] for fractionally charged colorless composite particles 
(fractons). 
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The U -hadron abundance in the interstellar gas strongly depends on the matter evolution 
in Galaxy, which is still not known to the extent, we need for our discussion. 

Indeed, in the opposite case of low density or of short time interval, when the condition 
(l32l is not valid, namely, at 



1 T tu 

n a < — ; — r~ = A-atom • T7 ■ — cm~ 3 < 10 3 for the case B, (35) 

fao{vv)t 300 K t ^ 



4 • 10 4 for the case A, 
10 3 for the case B , 
1.2 for the case C, 



where tu = 4 • 10 17 s is the age of the Universe, U -hadron abundance does not change its 
initial value. 

In principle, if in the course of evolution matter in the forming Galaxy was present 
during sufficiently long period (t ~ 10 9 yrs) within cold (T ~ 10 K) clouds with density 
n a ~ 10 3 cm -3 U -hadron abundance retains its primordial value for the cases A and B 
(/ = /), but falls down /„ = 5 • 10~ 9 in the case C, making this case close to the case B. 
The above argument may not imply all the U -hadrons to be initially present in cold clouds. 
They can pass through cold clouds and decrease their abundance in the case C at the stage of 
thermal instability, when cooling gas clouds, before they become gravitationally bound, are 
bound by the external pressure of the hot gas. Owing to their large inertia heavy [/-hadron 
atoms from the hot gas can penetrate much deeply inside the cloud and can be captured 
by it much more effectively, than ordinary atoms. Such mechanism can provide additional 
support for reduction of U -hadron abundance in the case C. The reduction of this abundance 
down to values, corresponding to the case A can be also provided by O-helium catalysis in 
the period of BBN. 

However, in particular, annihilation of [/-hadrons leads to multiple 7 production. If U- 
hadrons with relative abundance / annihilate at the redshift z, it should leave in the modern 
Universe a background 7 flux [15] 

F{E > E j = N^-f-n- s mod • c ^ 3 iq3/ (cm2 g ster) _^ 

of 7 quanta with energies E > E^ = 10 GeV/(l + z). The numerical values for 7 multi- 
plicity iV 7 are given in table 1 [15]. 



Energy fraction 


> 


> 0.1 GeV 


> 1 GeV 


> 10 GeV 


> 100 GeV 


A^ 7 


69 


62 


28 


2.4 


0.001 



Table 1: Multiplicities of 7 produced in the recombination of (QQ) pair with m = 250 
GeV for different energy intervals. 

So annihilation even as early as at z ~ 9 leads in the case C to the contribution into 
diffuse extragalactic gamma emission, exceeding the flux, measured by EGRET by three 
orders of magnitude. The latter can be approximated as 



F(E > E 7 ) ps 3- 10~ b I -^1 {cm 2 -s-ster)" 



where E = 451 MeV. 
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The above upper bound strongly restricts (/ < 1CT 9 ) the earliest abundance because 
of the consequent impossibility to reduce the primordial [/-hadron abundance by [/-hadron 
annihilation in low density objects. In the cases A and B annihilation in such objects should 
not take place, whereas annihilation within the dense objects, being opaque for 7 radiation, 
can avoid this constraint due to strong suppression of outgoing 7 flux. However, such con- 
straint nevertheless should arise for the period of dense objects' formation. For example, in 
the course of protostellar collapse hydrodynamical timescale tjj ~ 1/ yJnGp ~ 10 15 s/y^n 
exceeds the annihilation timescale [15] 

1 10 12 s /l/30\ 2 /_Z^\ 9 / 10 / m \ 11/1Q 

fn(av) fn \ a ) V300kJ 1,250 GeV7 

at n > 10 14 (i^J (i^) 4 (s^k) 975 (W&cv) 1175 ' wnere n is in cm ~ 3 - We consider 

homogeneous cloud with mass M has radius R 10 ^l/ c 3 m ^ > where M & = 2 • 10 33 

g is the Solar mass. It becomes opaque for 7 radiation, when this radius exceeds the mean 

free path / 7 ~ 10 26 cm/n at n > 3 • 10 10 cm -3 (xj^ ) 7 . As a result, for / as large as in 
the case C, rapid annihilation takes place when the collapsing matter is transparent for 7 
radiation and the EGRET constraint can not be avoided. The cases A and B are consistent 
with this constraint. 

_„ / M \2/9 

Note that at / < 5 • 10 \w^J > * e - ^ or a ^ tne considered cases energy release from 
U -hadron annihilation does not exceed the gravitational binding energy of the collapsing 
body. Therefore, U -hadron annihilation can not prevent the formation of dense objects but 
it can provide additional energy source, e.g. at early stages of evolution of first stars. Its 
burning is quite fast (few years) and its luminosity may be quite extreme, leading to a short 
inhibition of star formation [15]. Similar effects of dark matter annihilation were recently 
considered in [31]. 

Galactic blowing of CZ-baryon atoms polluting our Earth 
Since the condition Eq.(l35T) is valid for the disc interstellar gas, having the number density 
n g ~ 1 cm" 3 one can expect that the [7-hadron abundance in it can decrease relative to 
the primordial value only due to enrichment of this gas by the matter, which has passed 
through stars and had the suppressed U -hadron abundance according to Eq.(l33l). Taking the 
factor of such decrease of the order of the ratio of total masses of gas and stars in Galaxy 
f g ~ 10 -2 and accounting for the acceleration of the interstellar gas by Solar gravitational 
force, so that the infalling gas has velocity v g w 4.2 • 10 6 cm/s in vicinity of Earth's orbit, 
one obtains that the flux of [7-hadrons coming with interstellar gas should be of the order 
of [15] 

f fg n 9 v 9 _ 1 5 1Q -7_l r„„2 „ 
8vr ~ ' 10" 
where / is given by the Eq.(l34l. 

Presence of primordial {7-hadrons in the Universe should be reflected by their existence 
in Earth's atmosphere and ground. However, according to Eq.(l33l) (see discussion in Sec- 
tion 12.21 ) primordial terrestrial f7-hadron content should strongly decrease due to radiative 
recombination, so that the [/-hadron abundance in Earth is determined by the kinetic equi- 
librium between the incoming U -hadron flux and the rate of decrease of this abundance by 
different mechanisms. 



Iu = JJ9 2 9 ~ 1.5 • lO- 7 ^^ (sm 2 • s • ster) , (36) 
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In the successive analysis we'll concentrate our attention on the case, when the U 
baryon has charge +2, and ?7-hadrons are electrically neutral. In this case U baryons look 
like superheavy anomalous helium isotopes. 

Searches for anomalous helium were performed in series of experiments based on ac- 
celerator search [32], spectrometry technique [33] and laser spectroscopy [34]. From the 
experimental point of view an anomalous helium represents a favorable case, since it re- 
mains in the atmosphere whereas a normal helium is severely depleted in the terrestrial 
environment due to its light mass. 

The best upper limits on the anomalous helium were obtain in [34]. It was found by 
searching for a heavy helium isotope in the Earth's atmosphere that in the mass range 5 GeV 
- 10000 GeV the terrestrial abundance (the ratio of anomalous helium number to the total 
number of atoms in the Earth) of anomalous helium is less than (2—3) • 10~ 19 . The search in 
the atmosphere is reasonable because heavy gases are well mixed up to 80 km and because 
the heavy helium does not sink due to gravity deeply in the Earth and is homogeneously 
redistributed in the volume of the World Ocean at the timescale of 10 3 yr. 

The kinetic equations, describing evolution of anomalous helium and ^7-hadrons in 
matter have the form [15] 

dnjj 

= ju -n v -nfj- {av) - j gU (37) 

for J7-hadron number density n and 

dnrj 

-jf =Ju- n u- n u- {<rv) - j gU (38) 

for number density of anomalous helium njj- Here ju and j v take into account the income 
of, correspondingly, U -baryons and ?7-hadrons to considered region, the second terms on 
the right-hand-side of equations describe UU recombination and the terms j g u and j gU 
determine various mechanisms for outgoing fluxes, e.g. gravitationally driven sink of par- 
ticles. The latter effect is much stronger for [7-hadrons due to much lower mobility of 
U -baryon atoms. However, long range Coulomb like interaction prevents them from sink- 
ing, provided that its force exceeds the Earth's gravitational force. 

In order to compare these forces let's consider the World's Ocean as a thin shell of 
thickness L 4 • 10 5 cm with homogeneously distributed y charge, determined by distri- 
bution of [/-baryon atoms with concentration n. The y-field outside this shell according to 
Gauss' law is determined by 

2E y S = Aire y nSL, 

being equal to 

E y = 2ire y nL. 



In the result y force, exerting on [/-hadrons 

F y = e y E y , 

exceeds gravitational force for U -baryon atom concentration 



7 m 30 1 o 

n > 10 %nr v cm (39) 

250 GeV a y 
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Note that the mobility of [7-baryon atoms and U hadrons differs by 10 order of mag- 
nitude, what can lead to appearance of excessive y-charges within the limits of d39l >. One 
can expect that such excessive charges arise due to the effective slowing down of {7-baryon 
atoms in high altitude levels of Earth's atmosphere, which are transparent for U hadrons, as 
well as due to the 3 order of magnitude decrease of U hadrons when they enter the Earth's 
surface. 

Under the condition of neutrality, which is strongly protected by Coulomb-like y- 
interaction, all the corresponding parameters for J7-hadrons and [7-baryons in the Eqs.(l37T)- 
(I38T ) are equal, if Eq.(l39l> is valid. Provided that the timescale of mass exchange between 
the Ocean and atmosphere is much less than the timescale of sinking, sink terms can be 
neglected. 

The stationary solution of Eqs.(l37T)-(l38T) gives in this case 



n = J-±-, (40) 

V \ av ) 

where 

3u =3u=3~ 2 -^ = 1(rl2 IFIo cm_V1 (41) 

and (av) is given by the Eq.(l2TT). For j < 10~ 12 1Q t w cm~ 3 s~ 1 and (av) given by Eq. (|2TT) 
one obtains in water 

It corresponds to terrestrial [7-baryon abundance 

/«< io 2:1 ■ 1 



10 



-10 ' 



being below the above mentioned experimental upper limits for anomalous helium (f a < 
10~ 19 ) even for the case C with / = 2 • 10~ 6 . In air one has 



n^lO-y^cm- 3 . 

For example in a cubic room of 3m size there are nearly 27 thousand heavy hadrons. 

Note that putting formally in the Eq.(l40l the value of (av) given by the Eq.([T3l one 
obtains n < 6 • 10 2 cm~ 3 and f a < 6 • 10~ 20 , being still below the experimental upper 
limits for anomalous helium abundance. So the qualitative conclusion that recombination 
in dense matter can provide the sufficient decrease of this abundance avoiding the con- 
tradiction with the experimental constraints could be valid even in the absence of gauge 
y-charge and Coulomb-like y-field interaction for f7-hadrons. It looks like the hadronic 
recombination alone can be sufficiently effective in such decrease. However, if we take the 
value of (av) given by the Eq. (fT4l one obtains n by the factor of ^ ~ 5.5 larger and 
fa < 3.3 • 10~ 19 , what exceeds the experimental upper limits for anomalous helium abun- 
dance. Moreover, in the absence of y-attraction there is no dynamical mechanism, making 
the number densities of [/-baryons and [/-hadrons equal each other with high accuracy. So 
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Figure 1: Upper constraints on relative abundances of U-hadrons (/) following from the 
search for anomalous helium in the Earth, from EGRET, constraining possible annihilation 
effect on pregalactic stage, in comparison with prediction / in three cases (A,B,C). 

nothing seems to prevent in this case selecting and segregating U -baryons from £7-hadrons. 
Such segregation, being highly probable due to the large difference in the mobility of U- 
baryon atoms and U hadrons can lead to uncompensated excess of anomalous helium in the 
Earth, coming into contradiction with the experimental constraints. 

Similar result can be obtained for any planet, having atmosphere and Ocean, in which 
effective mass exchange between atmosphere and Ocean takes place. There is no such 
mass exchange in planets without atmosphere and Ocean (e.g. in Moon) and [7-hadron 
abundance in such planets is determined by the interplay of effects of incoming interstellar 
gas, UU recombination and slow sinking of U -hadrons to the centers of planets. 

Radiative recombination here considered is not able to save the case when free single 
positively charged U-hadrons (U ud) are present with noticeable fraction among U-hadronic 
relics. It is the very strong constraint on anomalous hydrogen (f a < 10~ 30 ) [35] what is 
hardly avoidable. So, the model with free stable +1 charged relics which are not bound in 
nuclear systems with charge > +2 can be discarded. 

To conclude this section we present the obtained constraints on the Fig.©. 

Cosmic rays experiments 
In addition to the possible traces of U-hadrons existence in the Earth, they can be manifested 
in cosmic rays. 

According to the arguments in previous section U -baryon abundance in the primary 
cosmic rays can be close to the primordial value /. It gives for case B 

/ = T ± ~ 3 • 1(T 9 . (42) 

n 

If [/-baryons have mostly the form (UUU), its fraction in cosmic ray helium component 
can reach in this case the value 
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Figure 2: Lines show the region of expected signal, their thickness reflects expected ac- 
curicy in experiment PAMELA. Upper thick curves relate to anomalous helium of 500 and 
750 GeV (a little upper and darker). The low thin curves relate to usual nuclei. 

which is accessible for cosmic rays experiments, such as RIM-PAMELA, being under run, 
and AMS 02 on International Space Station. 

Similar argument in the case C would give for this fraction ~ 2 • 10~ 5 , what may be 
already excluded by the existing data. However, it should be noted that the above estimation 
assumes significant contribution of particles from interstellar matter to cosmic rays. If 
cosmic ray particles are dominantly originated from the purely stellar matter, the decrease 
of [7-hadron abundance in stars would substantially reduce the primary [7-baryon fraction 
of cosmic rays. But in cosmic rays there are many different kind of particles. In order to 
differ U-hadrons from background events we can use the dependence of rigidity on velocity 
given on Figj2] 

The expected signal will strongly differ from background events in terms of relation be- 
tween velocity and rigidity (momentum related to the charge). In the experiment PAMELA 
velocity is measured with a good accuracy, what can lead to the picture Fig|2] 

Correlation between cosmic ray and large volume underground detectors' effects 
Inside large volume underground detectors (as Super Kamiokande) and in their vicinity U- 
hadron recombination should cause specific events ("spherical" energy release with zero 
total momentum or "wide cone" energy release with small total momentum), which could 
be clearly distinguished from the (energy release with high total momentum within "nar- 
row cone") effects of common atmospheric neutrino - nucleon-lepton chain (as well as of 
hypothetical WIMP annihilation in Sun and Earth) [15]. 

The absence of such events inside 22 kilotons of water in Super Kamiokande (SK) 
detector during 5 years of its operation would give the most severe constraint 

n < 10~ 3 cm -3 , 

corresponding to f a < 10~ 26 . For the considered type of anomalous helium such constraint 
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would be by 7 order of magnitude stronger, than the results of present direct searches and 3 
orders above our estimation in previous Section. 

However, this constraint assumes that distilled water in SK does still contain polluted 
heavy hadrons (as it may be untrue). Nevertheless even for pure water it may not be the case 
for the detector's container and its vicinity. The conservative limit follows from the condi- 
tion that the rate of [/-hadron recombination in the body of detector does not exceed the rate 
of processes, induced by atmospheric muons and neutrinos. The presence of clustered-like 
muons originated on the SK walls would be probably observed. 

High sensitivity of large volume detectors to the effects of U -hadron recombination 
together with the expected increase of volumes of such detectors up to 1 km 3 offer the pos- 
sibility of correlated search for cosmic ray U -hadrons and for effects of their recombination. 

During one year of operation a 1 km 3 detector could be sensitive to effects of recombi- 
nation at the [/-hadron number density n « 7 ■ 10~ 6 cm" 3 and f a « 7 ■ 10~ 29 , covering 
the whole possible range of these parameters, since this level of sensitivity corresponds 
to the residual concentration of primordial [/-hadrons, which can survive inside the Earth. 
The income of cosmic [/-hadrons and equilibrium between this income and recombination 
should lead to increase of effect, expected in large volume detectors. 

Even, if the income of anomalous helium with interstellar gas is completely suppressed, 
pollution of Earth by U -hadrons from primary cosmic rays is possible. The minimal effect 
of pollution by [/-hadron primary cosmic rays flux ly corresponds to the rate of increase 
of [/-hadron number density j ~ t^t> wnere Re ~ 6 • 10 8 cm is the Earth's radius. If 
incoming cosmic rays doubly charged [/-baryons after their slowing down in matter recom- 
bine with electrons we should take instead of Re the Ocean's thickness L « 4 • 10 5 cm that 
increases by 3 orders of magnitude the minimal flux and the minimal number of events, es- 
timated below. Equilibrium between this income rate and the rate of recombination should 
lead to N ~ jVt events of recombination inside the detector with volume V during its 
operation time t. 

For the minimal flux of cosmic ray U -hadrons, accessible to AMS 02 experiment during 
3 years of its operation (I m in ~ 10~ 9 / a ~ 4 • 10^ u 1(E), in the range of energy per 
nucleon 1 < E < 10 GeV) the minimal number of events expected in detector of volume 
V during time t is given by N m i n ~ 27r ^/" Vt. It gives about 3 events per 10 years in 
Super Kamiokande (V = 2.2 • 10 10 cm 3 ) and about 10 4 events in the 1 km 3 detector during 
one year of its operation. The noise of this rate is one order and half below the expected 
influence of atmospheric v^. 

The possibility of such correlation facilitates the search for anomalous helium in cosmic 
rays and for the effects of U -hadron recombination in the large volume detectors. 

The previous discussion assumed the lifetime of [/-quarks r exceeding the age of the 
Universe tjj- In the opposite case r < tjj all the primordial U hadrons should decay to the 
present time and the cosmic ray interaction may be the only source of cosmic and terrestrial 
U hadrons. 
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3 The case of a charge-asymmetry of U-quarks 



The model [15] admits that in the early Universe an antibaryon asymmetry for 4th genera- 
tion quarks can be generated [16,20,21]. Due to y-charge conservation U excess should be 
compensated by TV excess. We will focus our attention here to the case of y-charged quarks 
and neutrinos of 4th generation and follow [16] in our discussion. All the main results 
concerning observational effects, presented here, can be generalized for the case without 
y-interaction. 

[/-antibaryon density can be expressed through the modern dark matter density = 
k • Qcdm = 0.224 (k < 1), saturating it at A: = 1. It is convenient to relate the baryon 
(corresponding to f2& = 0.044) and U (TV) excess with the entropy density s, introducing 
Tb = fib/s and Tjj = njj/s = 3 • n^/s = 3 • r^. One obtains r& ~ 8 • 10 -11 and 
rjj, corresponding to U excess in the early Universe Kg = — rjj = 3 • (r^ — rjv) = 
10- 12 (350GeV/m[7) = lO" 12 /^, where S 5 = mt//350GeV. 

Primordial composite forms of 4th generation dark matter 
In the early Universe at temperatures highly above their masses U and TV were in thermo- 
dynamical equilibrium with relativistic plasma. It means that at T > vrijj (T > ttin) the 
excessive U (TV) were accompanied by UU (TVTV) pairs. 

Due to U excess frozen out concentration of deficit U -quarks is suppressed at T < mjj 
for k > 0.04 [21]. It decreases further exponentially first at T ~ l\j « a 2 Mjj /2 ~ 
35 5 GeV (where [15] a = C F a c = 4/3 • 0.144 w 0.19 and M v = m v /2 is the reduced 
mass), when the frozen out U quarks begin to bind with antiquarks U into charmonium-like 
state (UU) and annihilate. On this line U excess binds at T < Ijj by chromo-Coulomb 
forces dominantly into (UUU) anutium states with electric charge Za = —2 and mass 
m = 1.0555 TeV, while remaining free U anti-quarks and anti-diquarks (UU) form after 
QCD phase transition normal size hadrons (Uu) and (UUu). Then at T = Tqcd ~ 
150MeV additional suppression of remaining U -quark hadrons takes place in their hadronic 
collisions with £7-hadrons, in which (UU) states are formed and [/-quarks successively 
annihilate. 

Effect of TV excess in the suppression of deficit N takes place at T < m n f or k > 0.02 
[21]. At T ~ I NN = a 2 y M N /A ~ 15MeV (for a y = 1/30 and M N = 50GeV) due to 
y-interaction the frozen out N begin to bind with TV into charmonium-like states (TV TV) and 
annihilate. At T < I^u = o? y M^ )2 ~ 30MeV y-interaction causes binding of TV with 
U -hadrons (dominantly with anutium) but only at T ~ I^u /30 ~ lMeV this binding is 
not prevented by back reaction of y-photo-destruction. 

To the period of Standard Big Bang Nucleosynthesis (SBBN) U are dominantly bound 
in anutium with small fraction (~ 10~ 6 ) of neutral (Uu) and doubly charged (UUu) 
hadron states. The dominant fraction of anutium is bound by y-interaction with TV in 
(TVA~£~ ) "atomic" state. Owing to early decoupling of y-photons from relativistic plasma 
presence of y-radiation background does not influence SBBN processes [15, 16,20]. 

At T < I = Z 2 Z 2 He a 2 m He /2 w 1.6MeV the reaction A~^T + 4 He -► 7 + 
( 4 He ++ A~^) might take place, but it can go only after A He is formed in SBBN at 
T < lOOkeV and is effective only at T < T rHe ~ 7 /log (n^/n He ) w I Q /27 w 60keV, 
when the inverse reaction of photo-destruction cannot prevent it [14, 20, 22, 36]. In this 
period anutium is dominantly bound with TV. Since ru e = 0.1r& 3> r& = rfj/3, in this 
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reaction all free negatively charged particles are bound with helium [14, 20, 22, 36] and 
neutral Anti-Neutrino-O-helium (ANO-helium, ANOHe) ( 4 He ++ [NA~^}) "molecule" 
is produced with mass nioHe ~ m o ~ l^TeV. The size of this "molecule" is R ~ 
l/{Z/\ZH e otmHe) ~ 2 ■ 10~ 13 cm and it can play the role of a dark matter component and 
a nontrivial catalyzing role in nuclear transformations. 

In nuclear processes ANO-helium looks like an a particle with shielded electric charge. 
It can closely approach nuclei due to the absence of a Coulomb barrier and opens the way 
to form heavy nuclei in SBBN. This path of nuclear transformations involves the fraction 
of baryons not exceeding 1CT 7 [20] and it can not be excluded by observations. 

ANO-helium catalyzed processes 
As soon as ANO-helium is formed, it catalyzes annihilation of deficit JJ-hadrons and N. 
Charged [/-hadrons penetrate neutral ANO-helium, expel A He, bind with anutium and anni- 
hilate falling down the center of this bound system. The rate of this reaction is (av) = irR 2 , 
and an U excess k = 10~ 3 is sufficient to reduce the primordial abundance of (Uud) below 
the experimental upper limits. N capture rate is determined by the size of (NA) "atom" in 
ANO-helium and its annihilation is less effective. 

The size of ANO-helium is of the order of the size of A He and for a nucleus A with 
electric charge Z > 2 the size of the Bohr orbit for a (ZA) ion is less than the size of nu- 
cleus A. This means that while binding with a heavy nucleus A penetrates it and effectively 
interacts with a part of the nucleus with a size less than the corresponding Bohr orbit. This 
size corresponds to the size of 4 He, making O-helium the most bound (ZA)-atomic state. 

The cross section for A interaction with hadrons is suppressed by factor ~ {ph/PA) 2 ~ 
{r^/rh) 2 ps 10~ 4 /S 2 , where ph and pa are quark transverse momenta in normal hadrons 
and in anutium, respectively. Therefore anutium component of (ANOHe) can hardly be 
captured and bound with nucleus due to strong interaction. However, interaction of the 
4 He component of (ANOHe) with a nucleus can lead to a nuclear transformation due 
to the reaction + (AHe) Q + ^> provided that the masses of the initial and 

final nuclei satisfy the energy condition M(A, Z) + M(4, 2) - I Q > M(A + 4, Z + 2), 
where I Q = 1.6MeV is the binding energy of O-helium and M(4, 2) is the mass ofthe 4 tfe 
nucleus. The final nucleus is formed in the excited [a, M(A, Z)\ state, which can rapidly 
experience a- decay, giving rise to (ANOHe) regeneration and to effective quasi-elastic 
process of (AiVOi^e)-nucleus scattering. It leads to possible suppression of ANO-helium 
catalysis of nuclear transformations in matter. 

ANO-helium dark matter 
At T < T d ~ 1 keV energy and momentum transfer from baryons to ANO-helium 
njj(crv) (m p /m )t < 1 is not effective. Here a ps a Q ~ ttR 2 ps 10 _25 cm 2 . and 

v = ^2T/m p is baryon thermal velocity. Then ANO-helium gas decouples from plasma 
and radiation and plays the role of dark matter, which starts to dominate in the Universe at 
T RM = leV. 

The composite nature of ANO-helium makes it more close to warm dark matter. The 
total mass of (OHe) within the cosmological horizon in the period of decoupling is inde- 
pendent of 55 and given by 

Mod = -7^mpi(—-) Ps 2 • 10 g = lO y M . 

J-od J-od 
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O-helium is formed only at T Q = 60 keV and the total mass of OHe within cosmological 
horizon in the period of its creation is M a = M d(T /T d) 3 = 10 37 g. Though after 
decoupling Jeans mass in (OHe) gas falls down Mj ~ 3 • 10~ 14 M oc i one should expect 
strong suppression of fluctuations on scales M < M a as well as adiabatic damping of sound 
waves in RD plasma for scales M < M < M ^. It provides suppression of small scale 
structure in the considered model. This dark matter plays dominant role in formation of 
large scale structure at A; > 1/2. 

The first evident consequence of the proposed scenario is the inevitable presence of 
ANO-helium in terrestrial matter, which is opaque for (ANOHe) and stores all its in- 
falling flux. If its interaction with matter is dominantly quasi-elastic, this flux sinks down 
the center of Earth. If ANO-helium regeneration is not effective and A remains bound with 
heavy nucleus Z, anomalous isotope of Z — 2 element appears. This is the serious problem 
for the considered model. 

Even at k = 1 ANO-helium gives rise to less than 0.1 [20,21] of expected background 
events in XQC experiment [37], thus avoiding for all k < 1 severe constraints on Strongly 
Interacting Massive particles SIMPs obtained in [3] from the results of this experiment. 
In underground detectors (ANOHe) "molecules" are slowed down to thermal energies 
far below the threshold for direct dark matter detection. However, (ANOHe) destruction 
can result in observable effects. Therefore a special strategy in search for this form of 
dark matter is needed. An interesting possibility offers development of superfluid 3 He 
detector [38]. Due to high sensitivity to energy release above (Efh = 1 keV), operation of 
its actual few gram prototype can put severe constraints on a wide range of k and £5. 

At 10 -3 < k < 0.02 [7-baryon abundance is strongly suppressed [16,21], while the 
modest suppression of primordial N abundance does not exclude explanation of DAMA, 
HEAT and EGRET data in the framework of hypothesis of 4th neutrinos but makes the 
effect of ./V annihilation in Earth consistent with the experimental data. 

4 Signatures for U -hadrons in accelerator experiments 

Metastable [/-quark within a wide range of expected mass can be searched on LHC and 
Tavatron. In spite on that its mass can be quite close to that of t-quark, strategy of their 
search should be completely different, [/-quark in framework of considered model is 
metastable and will form metastable hadrons at accelerator contrary to i-quark. 

Detailed analysis of possibility of U-quark search requires quite deep understanding 
of physics of interaction between metastable U-hadrons and nucleons of matter. However, 
strategy of U-quark search can be described in general outline, by knowing mass spectrum 
of U-hadrons, (differential) cross sections of their production. LHC certainly will provide a 
better possibility for U-quark search than Tevatron. Cross section of U-quark production in 
pp-collisions at the center mass energy 14 TeV is presented on the Fig. [3] For comparison, 
cross sections of 4th generation leptons are shown too. Cross sections of U- and D- quarks 
does not virtually differ. 

Heavy metastable quarks will be produced with high transverse momentum pr, velocity 
less than speed of light. In general, simultaneous measurement of velocity and momentum 
enables to find the mass of particle. Information on ionization losses is, as a rule, not so 
good thereto. All these features are typical for any heavy particle, while there can be subtle 
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Figure 3: Cross sections of production of 4th generation particles (N, E, U, D) at LHC. 
Horizontal dashed line shows approximate level of sensitivity to be reached after first year 
of LHC operation. 

differences in the shapes of its angle-, py-distributions, defined by concrete model which it 
predicts. 

It is peculiarities of long-lived hadronic nature what can be of special importance for 
clean selection of events of U-quarks creation. U-quark can form a whole class of U- 
hadronic states which can be perceived as stable in condition of experiment contrary to their 
relics in Universe. However, as we pointed out, double, triple U-hadronic states cannot be 
virtually created in collider. Many other hadronic states whose lifetime is > 1CT 7 s should 
look like stable. In the Table 2 expected mass spectrum of U-hadrons, obtained with the 
help of code Pythia [39], is presented. 

The lower indexes in notation of U-hadrons in the Table 2 mean (iso)spin (/) of the light 
quark pair. From comparison of masses of different U-hadrons it follows that all I = 1 U- 
hadrons decay quickly emitting 7r-meson or 7-quantum, except (C/\ss)-state. In the right 
column the expected relative yields are present. Unstable 1 = 1 U-hadrons decay onto 
respective 1 = states, increasing their yields. 

Firstly one makes a few notes. There are two mesonic states being quasi-degenerated 
in mass: (Uu) and (Ud) (we skip here discussion of strange U-hadrons). In interaction 
with medium composed of u and d quarks transformations of U-hadrons into those ones 
containing u and d are preferable (as it is the case in early Universe). From these it follows, 
that created pair of U U quarks will fly out from the vertex of pp-collision in form of U- 
hadron with positive charge in about 60% of such events and with neutral charge in 40% and 
in form of anti-U-hadron with negative charge in 60% and neutral in 40%. After traveling 
through detectors a few nuclear lengths from vertex, U-hadron will transform in (roughly) 
100% to positively charged hadron (Uud) whereas anti-U-hadron will transform in 50% to 
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Table 2: Mass spectrum and relative yields in LHC for U-hadrons. The same is for charged 
conjugated states. 





Difference between the 
masses of U-hadron 
and U-quark, GeV 


Expected yields (in the right 
columns the yields of long-lived 
states are given) 


{Uuf ,{Ud} + 


0.330 


39.5(3)%, 39.7(3)% 


{usy 


0.500 


11.6(2)% 


{Uud} + 


0.579 


5.3(1)% 


7.7(1)% 


{Uuu}+ + ,{Uud}t , 
{Udd}^ 


0.771 


0.76(4)%, 
0.86(5)%, 
0.79(4)% 




{Usu} + ,{Usd} u 


0.805 


0.65(4)%, 
0.65(4)% 


1.51(6)% 


{Uau}f , {Usdj'i 


0.930 


0.09(2)%, 
0.12(2)% 




{u ss y (1) 


1.098 


0.005(4)% 



negatively charged U-hadron (Ud) and in 50% to neutral U-hadron (U u). 

This feature will enable to discriminate the considered model of U-quarks from variety 
of alternative models, predicting new heavy stable particles. 

Note that if the mass of Higgs boson exceeds 2m, its decay channel into the pair of 
stable QQ will dominate over the ti, 2W, 2Z and invisible channel to neutrino pair of 4th 
generation [40]. It may be important for the strategy of heavy Higgs searches. 

5 Conclusion 

To conclude, the existence of hidden stable or metastable quark of 4th generation can be 
compatible with the severe experimental constraints on the abundance of anomalous iso- 
topes in Earths atmosphere and ground and in cosmic rays, even if the lifetime of such 
quark exceeds the age of the Universe. Though the primordial abundance / = r^jr^ of 
hadrons, containing such quark (and antiquark) can be hardly less than / ~ 10~ 10 in case 
of charge symmetry, their primordial content can strongly decrease in dense astrophysi- 
cal objects (in the Earth, in particular) owing to the process of recombination, in which 
hadron, containing quark, and hadron, containing antiquark, produce unstable charmonium- 
like quark-antiquark state. 

To make such decrease effective, the equal number density of quark- and antiquark- 
containing hadrons should be preserved. It appeals to a dynamical mechanism, preventing 
segregation of quark- and antiquark- containing hadrons. Such mechanism, simultaneously 
providing strict charge symmetry of quarks and antiquarks, naturally arises, if the 4th gen- 
eration possesses new strictly conserved U(l) gauge (y-) charge. Coulomb-like y-charge 
long range force between quarks and antiquarks naturally preserves equal number densities 
for corresponding hadrons and dynamically supports the condition of y-charge neutrality. 

It was shown in the present paper that if [/-quark is the lightest quark of the 4th genera- 
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tion, and the lightest free [/-hadrons are doubly charged (UUU)- and (L r f/u)-baryons and 
electrically neutral (Uu) -meson, the predicted abundance of anomalous helium in Earths 
atmosphere and ground as well as in cosmic rays is below the existing experimental con- 
straints but can be within the reach for the experimental search in future. To realize this 
possibility nuclear binding of all the (Uud) -baryons with primordial helium is needed, con- 
verting potentially dangerous form of anomalous hydrogen into less dangerous anomalous 
lithium. Then the whole cosmic astrophysics and present history of these relics are puzzling 
and surprising, but nearly escaping all present bounds. 

Searches for anomalous isotopes in cosmic rays and at accelerators were performed 
during last years. Stable doubly charged U baryons offer challenge for cosmic ray and 
accelerator experimental search as well as for increase of sensitivity in searches for anoma- 
lous helium. In particular, they seem to be of evident interest for cosmic ray experiments, 
such as PAMELA and AMS02. +2 charged U baryons represent the low Z/A anomalous 
helium component of cosmic rays, whereas —2 charged U baryons look like anomalous 
antihelium nuclei. In the baryon asymmetrical Universe the predicted amount of primor- 
dial single (Uud) baryons is exponentially small, whereas their secondary fluxes originated 
from cosmic ray interaction with the galactic matter are predicted at the level, few order of 
magnitude below the expected sensitivity of future cosmic ray experiments. The same is 
true for cosmic ray +2 charged U baryons, if U -quark lifetime is less than the age of the 
Universe and primordial U baryons do not survive to the present time. 

The models of quark interactions favor isoscalar (Uud) baryon to be the lightest among 
the 4th generation baryons (provided that U quark is lighter, than D quark, what also may 
not be the case). If the lightest U -hadrons have electric charge +1 and survive to the present 
time, their abundance in Earth would exceed the experimental constraint on anomalous 
hydrogen. This may be rather general case for the lightest hadrons of the 4th generation. 
To avoid this problem of anomalous hydrogen overproduction the lightest quark of the 4th 
generation should have the lifetime, less than the age of the Universe. Another possible 
solution of this problem, using double and triple U baryons (UUU) and catalysis of (Uud) 
annihilation in atom-like bound systems 4 He(UUU) is considered in [15]. 

However short-living are these quarks on the cosmological timescale in a very wide 
range of lifetimes they should behave as stable in accelerator experiments. For example, 
with an intermediate scale of about 10 11 GeV (as in supersymmetry models [41]) the ex- 
pected lifetime of U-(or D-) quark ~ 10 6 years is much less than the age of the Universe 
but such quark is practically stable in any collider experiments. 

First year operation of the accelerator LHC has good discovery potential for U(D)- 
quarks with mass up to 1.5 TeV. U-hadrons born at accelerator will distinguish oneself by 
high p t , low velocity, by effect of a charge flipping during their propagation through the 
detectors. All these features enable strongly to increase efficiency of event selection from 
not only background but also from alternative hypothesis. 

In the present work we studied effects of 4th generation having restricted our analysis by 
the processes with 4th generation quarks and antiquarks. However, as we have mentioned 
in the Introduction in the considered approach absolutely stable neutrino of 4th generation 
with mass about 50 GeV also bears y-charge. The selfconsistent treatment of the cosmo- 
logical evolution and astrophysical effects of y-charge plasma of neutrinos, antineutrinos, 
quarks and antiquarks of 4th generation in charge symmetric case will be the subject of 
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special studies. An attempt of such a treatment has been undertaken in the case of charge 
asymmetry, described in this paper. 

We believe that a tiny trace of heavy hadrons as anomalous helium and stable neutral 
O— Helium and meson£j]may be hidden at a low level in our Universe ~ 1CP 10 — 10~ 9 ) 
and even at much lower level here in our terrestrial matter a density ^ ~ 10~ 23 in case of 
charge symmetry. There are good reasons to bound the 4th quark mass below TeV energy. 
Therefore the mass window and relic density is quite narrow and well defined, open to a 
final test. 

In case of charge asymmetry of 4th generation quarks, a nontrivial solution of the 
problem of dark matter (DM) can be provided due to neutral O— Helium-like U-hadrons 
states (ANO-helium in case of y-interaction existence). Such candidates to DM have 
many interesting implications in BBN, large scale structure of Universe and physics of 
DM [16,20-26]. It should catalyze new types of nuclear transformations, reminding al- 
chemists' dream on the philosopher's stone. It challenges direct search for species of such 
composite dark matter and its constituents. A very low probability for their existence is 
strongly compensated by the expectation value of their discovery. 
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